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We impose new severe constraints on the self-interactions of fermionic asymmetric dark matter 
based on observations of nearby old neutron stars. WIMP self-interactions mediated by Yukawa- 
type interactions can lower significantly the number of WIMPs necessary for gravitational collapse 
of the WIMP population accumulated in a neutron star. Even nearby neutron stars located at 
regions of low dark matter density can accrete sufficient number of WIMPs that can potentially 
collapse, form a mini black hole, and destroy the host star. Based on this, we derive constraints on 
the WIMP self-interactions which in some cases are by several orders of magnitude stricter than the 
ones from the bullet cluster (which are currently considered the most stringent). 
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It is highly possible that dark matter (DM) is in the 
form of Weakly Interacting Massive Particles (WIMPs) . 
Although the cross section in WIMP-nucleon interac- 
tions is severely constrained in underground detectors, we 
have little information regarding self-interactions among 
WIMPs other than gravitational attraction. The possi- 
bility of WIMP self-interactions has been studied in the 
past pQ, and in principle there is no argument against 
it. On the contrary WIMP self-interactions have been 
proposed as a solution of the discrepancy between sim- 
ulations of collisionless DM and observations of galax- 
ies [2J [3J. Simulations of collisionless cold DM predict 
a cusped profile with a steep power law-like behavior in 
the DM density [3]. This picture is contradicted by ob- 
servations of the rotation curves of gas-rich dwarf galax- 
ies. Such galaxies are dominated by DM and therefore 
it is expected to give more reliable information regard- 
ing the DM density profile. Observations of such dwarf 
galaxies suggest that their DM halos have an approxi- 
mately constant density rather than a power law scal- 
ing. Generally speaking self-interactions can alleviate 
this core/cusp problem. Collisions among WIMPs can 
potentially flatten out the cusp, thus creating the almost 
constant density favored by observations. Such WIMP 
self-interactions should in general be strong enough in or- 
der to have frequent WIMP- WIMP collisions at the core, 
but they have to be short range interactions in order to 
avoid changing large scale structure [5]. In addition, for 
such a scenario to be realized, WIMP self-annihilations 
should be either absent or restricted. 

Although asymmetric dark mater is not a new idea [BJ- 
9 , it has recently attracted a lot of interest p.0 27 . This 
is because it can be an alternative to thermally produced 
WIMP DM, it can link the asymmetry in the dark sector 
to the baryonic asymmetry, and it can possibly reconcile 
discrepancies among underground direct search experi- 
ments through interference of different channels [25J • 
In an asymmetric DM scenario, annihilations are rare due 
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to the fact that most of anti- WIMPs have already been 
annihilated in the past. Therefore self- interacting asym- 
metric DM is quite an attractive possibility for resolving 
the cusp/core problem. 

In this letter we impose model-independent con- 
straints on self-interacting asymmetric DM candidates 
with Yukawa-type interactions based on neutron star ob- 
servations. Compact stars such as neutron stars and 
white dwarfs can potentially constrain WIMP masses 
and cross sections with nuclei [130441). For asymmet- 
ric DM, such constraints are based on the fact that if 
enough WIMPs have been accreted and trapped within 
the compact star, they can collapse forming a black hole 
that eventually destroy the host star. Therefore the ex- 
istence of relatively old compact stars in our Galaxy in 
regions of sufficiently rich DM density can exclude a lot 
of phase space of asymmetric DM candidates. Based 
on such arguments, spin-dependent interactions of asym- 
metric WIMPs [35], and spin-(in)dependent interactions 
of bosonic asymmetric DM [?D1 2T] have been severely 
constrained. However, in the case of fermionic DM, con- 
straints based on the above argument are very weak. In 
order for WIMPs to collapse and overcome the WIMP 
Fermi pressure, a large number of trapped WIMPs is re- 
quired, making the constraints valid only for very large 
(and uninteresting) WIMP masses. However, the exis- 
tence of attractive self-interactions tends to alleviate the 
effect of Fermi pressure, changing drastically the number 
of WIMPs required for collapse. In this letter we assume 
that DM WIMPs (of mass m) self-interact via a Yukawa- 
type interaction of the form a4>ipip- If 4> (with a mass 
/i) is a scalar, such self-interactions are always attrac- 
tive, inducing a short range "dark force" with a potential 
V(r) — — aexp[— [ir}/r. This is the type of WIMP-self- 
interactions that have been proposed as an explanation of 
the flatness of DM density in dwarf galaxies [3] , or in or- 
der to fit the direct DM search experimental results [42j . 

WIMPs passing through a neutron star can scatter and 
potentially be trapped inside. The accretion of WIMPs 
onto a typical 1.4Mq 10 km neutron star taking into 
account relativistic effects has been calculated in [3TJ [35] . 
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The number of accumulated WIMPs is 

(i) 

where v is the WIMP average velocity, tg is the time in 
billion years, and / is a factor which is 1 if the WIMP- 
nucleon cross section a > 10 cm , and cr/(lCT 45 cni 2 ) 
if a < 10~ 45 cm 2 . Once the WIMPs get trapped, they 
thermalize with nucleons acquiring the same temperature 
in a very short time interval [35] and follow a Fermi- 
Dirac distribution. The WIMP thermal radius can be 
determined by applying the virial theorem for a random 
WIMP i 

o/^n 8 „ 2 GNm 2 e-" r « _ ar .\ 
2(E k ) = -irGpmr 2 + + ( \ a + ape ^ ) , 
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where (Ek) is the average kinetic energy of a single 
WIMP, N is the total number of WIMPs, p is the core 
mass density of the neutron star, r is the WIMP ther- 
mal radius, and r,j is the distance between WIMP i 
and WIMP j. The sum runs over all WIMPs (except 
i). All quantities have been averaged over time and i. 
The virial theorem for a generic potential of the form 
V(r) = r n e- arP is 2{E k ) = n(V) - ap( r PV) [43 . It can 
be easily seen that the first two terms on the right hand 
side correspond to the gravitational potential energy due 
to the core of the neutron star and the self-gravitation 
of WIMPs respectively. The last two terms correspond 
to the average attractive Yukawa potential energy. Since 
most of the WIMPs are concentrated within the thermal 
radius r, if x is the mean distance between two closest 
neighbor WIMPs, x = r/N 1 / 3 . At the beginning, as 
WIMPs start accumulating in the star, the first term of 
Eq. ([2| is the dominant one. As N increases, the other 
terms start to become more important. It is clear that 
Yukawa forces can play a significant role if px << 1 or 

pr/N 1/3 « 1. (3) 

In this case, the exponential will be roughly close to ~ 1 
and the Yukawa force between two WIMPs will scale 
as a Coulomb one. In order to determine the thermal 
radius, one has to solve Eq. (pi for r keeping in mind 
that the kinetic energy Ek = (3/2)fcsT for nondegener- 
ate WIMPs, or E k = (9tt/4) 2 / 3 7V 2 / 3 /(2mr 2 ) (for non- 
relativistic degenerate WIMPs) once nA 3 > 2, where n 
is the WIMP number density, and A = h/ \f2ixmkBT . 
We assumed for simplicity that WIMPs have two degen- 
erate states (spin- up spin-down). The generalization to 
a higher number of states is straightforward. For typi- 
cal parameters of a neutron star (and a thermal radius 
r th = [9/c B T/(87rGpm)] 1 /2 = 220cm/ v /m 40J), WIMPs 
become degenerate once their number is N^ eg = 6 x 10 35 . 
There are two different possibilities as the number of 
trapped WIMPs increases depending on the parameters 
of the Yukawa forces. Yukawa self-interactions might 
dominate Eq. ^ either before or after the degenerate 
limit has been reached. 



Self- attraction after degeneracy. The kinetic energy is 
the one of a nonrelativistic degenerate gas, and therefore 
the thermal radius is 

(4) 

We have used a typical value p = 5 x 10 38 GeV/cm 3 . 
There are two possibilities here: Either self-attraction 
(3rd and 4th terms in Eq. (|2|) dominates after Eq. ^ 
is satisfied, or not. In the first case, a comparison be- 
tween the 1st and the sum of the 3rd and 4th term of 
Eq. ([2]) shows that self-attraction takes place after the 
saturation of Eq. ^ condition, if a < 4.38/i/m. At this 
point we define the quantity y = /ir/iV 1 / 3 which will be 
more convenient to work with instead of dealing with r. 
The particles become relativistic once the momentum be- 
comes comparable to the mass {Q-K/Af^N 1 / 3 ^ ~ m, or 
Vrd = (97r/4) 1 / 3 ^i/m. Now we have to estimate the po- 
tential energy due to the Yukawa interactions for y < 1. 
The energy is 

3 J 

where i denotes a specific particle in question and j runs 
over all the rest. Apparently if y > 1, due to the exponen- 
tial decay of the potential with separation distance be- 
tween particles, the most significant contribution comes 
from the nearest neighbor. However if y < 1, not only 
there is no exponential suppression for the nearest neigh- 
bor, but also more distant particles can contribute as 
long as rij < It is easy to show that for y < 1, 

the number of particles inside a sphere with radius 1/p 
is 1/y 3 . This is the number of particles that will con- 
tribute significantly to the potential energy. If we take 
an approximate average (I/ty,-) ~ p,, the Yukawa poten- 
tial energy is Ei = ap/y 3 . 

As we mentioned above, for a < 4.38/i/m, self- 
attraction dominates after Eq. ^ has been saturated. 
There are two possible cases here. The first one is to 
have self-attraction dominating before the particles be- 
come relativistic. In this case the Yukawa self-attraction 
overcomes nonrelativistic Fermi pressure once 

ap > (9 1 \ V3 jl_ (6) 
y 3 \ 4 / my 2 ' 

This happens when y < y\ — am/3.65p. As the WIMPs 
keep falling, they become eventually relativistic. In order 
to further collapse, self-attraction should win over the 
relativistic Fermi pressure. This means 

ap /9vr\ 1/3 p p 2 

—q— > ( —r ) <^a>7— . (7) 

Vrel V 4 J yrel m 2 

Therefore for 7p 2 /m 2 < a < 4.38/i/m, it is guaranteed 
that WIMPs can gravitationally collapse forming a black 



3 




FIG. 1: Excluded regions (in light blue) for the WIMP mass 
m and the mediator mass /i in logarithmic scales of GeV for 
three different Yukawa couplings a = 10~ 5 , 0.01, and 0.1 
from top to bottom. Light brown regions are also excluded 
by constraints on WIMP self interactions based on the bullet 
cluster. The larger exclusion region in plot 1 and the "tri- 
angle" regions in plots 2 and 3 correspond to gravitational 
collapse initiated before WIMPs have become degenerate. 



hole so long y < y\. This condition on y can be trans- 
formed to a condition on the number of accreted WIMPs 
through Eq. (4). In this letter we are going to consider 
limits derived only from nearby old neutron stars with 
DM densities of pdm = 0.3GeV/cm 3 . Globular clusters 
can give stricter constraints, but we have chosen here to 
present results that do not depend on the not well known 



DM density in globular clusters, but rather on the unam- 
biguous low densities of nearby neutron stars. There is 
one extra constraint that we have to consider. If the 
above conditions are satisfied, WIMPs will form a black 
hole. However, as it was estimated in [3D], only black 
holes with masses larger than 5.7 x 10 36 GeV can even- 
tually destroy the host star. Smaller black holes decay 
due to the fact that Hawking evaporation dominates over 
the black hole accretion of matter in the neutron star. 

The case where Yukawa forces dominate at a y < y re i 
can also in principle lead to formation of a black hole. 
However, we found that this happens for a < 7/i 2 /m 2 . 
Practically for such small couplings, there is not enough 
phase space, meaning that accretion of WIMPs in nearby 
neutron stars in billions years is not enough to meet the 
condition for collapse. 

We move now to the case where a > 4.38/i/m. This 
means that Yukawa self-attraction becomes important 
before y = 1. This is distinctively different from the 
previous case because the Yukawa potential here is still 
exponentially suppressed (and has 1/y instead of 1/y 3 
dependence). The virial theorem Eq. ([2| dictates the 
thermal radius of the WIMP sphere. Ignoring WIMP 
self- gravitation which is small for our case of interest, 
Eq. Q can be rewritten as 



2/3 



m y 2 



a 1 



8 ^1 IK „T?/3 



y 2 ■ 



(8) 

Note that in this case we deal with y > 1 and therefore 
we have approximated the sum over all WIMPs in Eq. ^ 
simply by the closest neighbor. There are up to three so- 
lutions for this equation in the case where a > 4.38/i/m. 
One of them lies at y < 1 and the rest at y > 1. As N 
increases, the two solutions with y > 1 come closer to 
each other. Eventually at a critical value N cr i t , the two 
solutions coincide and then disappear leaving only the 
solution with y < 1. We have solved numerically Eq. ^ 
to identify N cr u as well as the value of y that is a solution 
when N is slighly larger than N cr i t . We have compared 
the Yukawa potential energy to the Fermi one, and as 
before the constraint that we derive is that a > 7/i 2 /ra 2 . 
We should also mention here that the same result can be 
obtained by minimizing the energy of the system. One 
can easily check that the minimization of the energy of 
the system gives the virial Eq. ([8|. 

In all of the above constraints, we have implicitly 
assume that y re i < 1 (or equivalently [i < 0.52m) 
i.e. WIMPs become relativistic when already Yukawa 
forces are not exponentially suppressed. In the opposite 
case where y re i > 1, WIMPs become relativistic before 
Yukawa forces become effective. In this case we find that 
there is a small phase space to be excluded, but it is 
rather uninteresting since it lies on a region excluded al- 
ready by constraints on WIMP self-interactions based on 
the bullet cluster. This case can also give constraints on 
strongly coupled interactions (a > 1) but in this letter 
we focus on perturbative couplings a < 1. 
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Self- attraction before degeneracy. Depending on the 
parameters, it is possible for Yukawa self-attractions to 
become important before the WIMPs have reached the 
degeneracy limit which corresponds to y^g — 1.3 x 
10 4 /j,/y / m (both to and fi measured in GeV). This hap- 
pens when the last two terms of Eq. ^ become com- 
parable to the first one. We have identified numerically 
the parameter region that this happens and we have es- 
timated the number of WIMPs N s needed for that. Note 
that this number has to be smaller than Nd eg - In such a 
case, it is possible to form a WIMP black hole if y^ eg < 1 
and the Yukawa potential dominates over nonrelativistic 
and relativistic Fermi pressure. This leads to the follow- 
ing constraints respectively 

2 2 

a> (9tt/4) 2 / 3 ^y deg ~4,8 x 10 4 ^ & a > 7^, 

(9) 

where for all cases of interest the first one is much tighter 
and therefore sufficient for collapse. Once more in order 
to impose actual constraints the WIMP black hole must 
be larger than 5.7x 10 36 GeV in order to grow and destroy 
the host star. 

In Fig. 1 we present all the exclusion regions based on 
the existence of old neutron stars in the vicinity of the 
Earth where the DM density is p = 0.3 GeV/cm 3 , such as 
J0437-4715 05 and J0108-1431 05] (140 pc and 130 pc 



from the Earth, respectively). The fact that these stars 
have not been destroyed by a black hole after accreting 
WIMPs for at least one billion years, can impose severe 
constraints on the possible WIMP mass and the mass of 
the mediator that is exchanged. The constraints are valid 
as long as the WIMP-nucleon cross section is not smaller 
than the a cr u — 10 _45 cm 2 . However, we found that 
the exclusion regions do not change significantly with 
a as low as 10 _48 cm 2 . Old neutron stars in globular 
clusters can potentially give stronger constraints, albeit 
the arbitrariness of the local DM density. In addition, 
neutron stars in binary systems might increase by a fac- 
tor of few the dark matter accretion 08j. In Fig. 1 we 
present new constraints not excluded by the bullet cluster 
observations. Using the WIMP- WIMP cross section for 
Yukawa interactions [3] with a relative velocity of 4500 
km/s relevant for the bullet cluster, we exclude a whole 
new region for self-interactions (for asymmetric WIMPs) . 
Comparatively, if one were to estimate how stricter are 
the constraints derived here, a WIMP of mass 1 TeV, 
exchanging a scalar of 100 MeV with a — 0.01, would 
lead to a WIMP- WIMP cross section exclusion of 10 
orders of magnitude lower than the bullet cluster limit 
o S eif/m < 2 x 10^ 24 cm 2 /GeV 03] . 

I would like to thank P. Tinyakov for useful comments 
on the manuscript. 
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